Vibrationally-mediated molecular transistors 
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We investigate the steady-state electronic transport through a suspended dimer molecule coupled 
to leads. When strongly coupled to a vibrational mode, the electron transport is enhanced at the 
phonon resonant frequency and higher-order resonances. The temperature and bias determines 
the nature of the phonon-assisted resonances, with clear absorption and emission peaks. The strong 
coupling also induces a Frank-Condon-like blockade, suppressing the current between the resonances. 
We compare an analytical polaron transformation method to two exact numerical methods: the 
Hierarchy equations of motion and an exact diagonalization in the Fock basis. In the steady- 
state, our two numerical results are an exact match and qualitatively reflect the main features of 
the polaron treatment. Our results indicate the possibility of a new type of molecular transistor or 
sensor where the current can be extremely sensitive to small changes in the energies of the electronic 
states in the dimer. 
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I. INTRODUCTION 

Electron transport through real, and artificial, sus- 
pended molecular systems has revealed a rich tapestry 
of physical effects 1 - - — . Nano-mechanical systems of- 
fer the promise of observing quantum effects in mas- 
sive objects^—, which can be cooled, driven, measured 
and manipulated using electron transport 7 -* 1 ^ - — . Sin- 
gle molecule electronics offers similar promises, with the 
benefit of stronger coupling of vibrational and oscilla- 
tory modes to the transport process, and a large range 
of practical applications 2 ^— . Current efforts on single- 
molecule electronics include the exploration of the elec- 
tronic level structure and its effect on electronic transport 
together with the development of electronic devices and 
applications^ 5 .. 

Since the performance of molecular electronic devices 
depends on electron transfer between molecules, it ulti- 
mately requires a good understanding of coherent trans- 
port and dynamics. Molecular electronics has a variety of 
predecessors in chemistry and chemical reactions, known 
as electron-transfer reaction a 36 ' 37 . One particular exam- 
ple, the donor-bridge-acceptor (DBA) system has been 
studied for a long time and used to probe the mechanisms 
of charge transfer 38 . In this system, the donor molecule 
donates an electron to the acceptor via the third non- 
rigid or rigid molecule called a "bridge" (i.e., redox reac- 
tion). In the lab it is possible to construct a molecular 
electronics device by attaching electrodes to the donor 
and acceptor molecules of a mixed- valence dimer and ap- 
ply a bias voltage to cause electron transport. Aviram 
and Ratner~ made the first metal-DBA-metal system 
in analogy to p-n junctions. The transfer rate of inco- 
herent transfer can be explained within the well-known 
Forster and Dexter theory 4 ^—, and until recently, co- 
herent transport had been regarded to play only a small 
role in charge-transfer dynamics of molecules. However, 
this view has been changed in the last decade and some 
experiments indicate that coherent transport is indeed 



importan t 44 ' 45 . 

In a slightly different context, recent experiments in 
photosynthetic complexes have shown that excitonic co- 
herence lasts surprisingly long at non-zero temperature: 
up to 300 fs at room temperature 4 ^— _ ft nas recently 
been propose d 49 ' 50 that the vibrational modes in the pro- 
tein environment may also play a role in the observed 
oscillations, long coherence time, and high transport ef- 
ficiency. Therefore, an understanding of the intricacies 
of both coherent electron transport and coherent exciton 
transport transport in molecules, coupled to vibrational 
modes, is vital22i£i£ir£Z. 

In this paper, we examine a model of electron 
transport through a two-site system coupled linearly 
to a single vibrational mode. Motivated by recent 
progress in single-molecule experiments, we study the 
case when the electron-phonon coupling can be strong 
(non-perturbative), and where the phonons can have a 
life-time much longer than the electronic transport pro- 
cess. Throughout we use the language of a single mixed- 
valence dimer molecule coupled to electronic leads. One 
electrode (source) is attached to one site of a dimer 
molecule (donor) and another electrode (drain) is at- 
tached to the other site (acceptor). This system is analo- 
gous to the DBA system, but instead of incoherent trans- 
port, we include the possibility of coherent electron trans- 
port and analyze steady-state dynamics and current noise 
power spectra. 

We begin with a semi-analytical treatment based 
on a polaron-transform, analogous to a single- 
mode version of the traditional "non-interacting blip 
approximation" (NIBA) 58 ' 59 , which can give us some in- 
sight on a regime not easily accessible with numerical 
methods (non-equilibrium transport, un-damped vibra- 
tional mode and strong coupling) . Such an approach has 
been used elsewhere to study the transport through sin- 
gle and double quantum dot o 13 ' 14 ' 60 ' 61 , photosynthetic 
complexes 62 (and similar, though slightly different, tech- 
niques were used i n 49 i 50 ' 63 ), and molecules. Our deriva- 
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tion follows in the same vein as these earlier works. We 
follow this with two exact numerical treatments, one 
based on the Hierarchical equations of motio n 64 ! 65 with 
a Lorentz bath, the other an exact diagonalization of the 
Hamiltonian in the Fock basis. 

Using these various methods we show that in certain 
conditions the vibrational mode coupling can strongly 
enhance electron transport in a selective manner, sug- 
gesting the possibility of using vibrationally-assisted elec- 
tronic transitions as a highly sensitive transistor. 



II. HAMILTONIAN 

Consider a dimer molecule (denoted below by sites 
left (L) and right (R)) with each site coupled to a lead 
(fermionic reservoir) so that the dimer functions as a 
bridge. A vibrational mode, e.g., a flexural or dila- 
tional mode of the molecule (or bridge in the D-B-A sys- 
tem), is coupled to the electrons in the left and right 
sites of the dimer. The idea of strong coupling to one 
vibrational mode is not far-fetche d 33 i 61 : Such coupling 
has been shown experimentally in some nano-electronics 
devices 6 ^ -. 

The Hamiltonian of this system, which is analogous to 
the Hubbard-Peierls model in the large on-site repulsion 
limit and with spin suppressed, can be written as 



H — Hq + Ht + H\ 



where 



Hq — -ffmol H 
Ht = -fftunn 



H 



couple 



H v = VLa\ds + V R d} D a R + V L d\a L + V R a R d D 



H 



1 t if 

mol = 2 £ i fl L a i + rQR 

H osc = hhJmb^b 
#icads = ^2 {j^kdl^d^k + fruj? d} D k d D + 

k 

#tunn = T c (a\a R + a R ai)j 
-ffcoupio = 9l (6 f + b) a\a L - g R (V + b) a R a R . 



(1) 
(2) 

(3) 
(4) 
(5) 

(6) 
(7) 



Here Hq is the free Hamiltonian of the dimer H mo u the 
oscillator H osc , and the source and drain reservoirs 7?i oa ds 
attached to the left (right) sites, with resonant frequen- 



cies U) L , UJ R , OJr, 



J k > 



and 



J k ' 



respectively. Also, T c is 



the transmission coefficient, and gL(R) is the electron- 
vibrational-mode coupling coefficient. The annihilation 
(creation) operators for the electron in the molecules are 

a L(R) ( a l(_R))i f° r the lead reservoirs are cis(d) ( a s(£>))' 
and for the oscillator is b (&*). The interaction Hamil- 
tonian has three parts: lead electrode-molecule coupling 
Hy , electron tunneling between the two sites in the dimer 



molecule H tunn , and the linear coupling between the 
dimer and the vibrational mode iJ C oupio- The coupling 
between the electrons and the vibrational mode effec- 
tively modify the frequency of the electrons in the left 
and right sites of the molecule as ujl(r.) ~^ w l(_r)' wnere 
= u l(r) + 9l(R) \b' + b), while the tunneling bar- 
rier is unaffected by the coupling. 

Here we adopt the notation, for brevity, 



t t 



a R a L 



(8) 



In the limit of large Coulomb repulsion to double occu- 
pation (a term not explicitly included in the Hamiltonian 
for simplicity), only a single electron occupies the entire 
molecule dimer at once. For the moment, we retain the 
notation for the general case, but will soon explicitly use 
this assumption. 

The above Hamiltonian is the most general and exact 
form. To find quasi-analytical results, this Hamiltonian 
can be transformed using the polaron transformation to 
eliminate the interaction term, but which in turn modifies 
the tunneling term Ht- 



H t - T c (st^t + sX) 



where 



X = exp 



A f = exp -z (b - P 



(9) 

(10) 
(11) 



with z = g/cj m , g = Ql = gR- Note that we have ignored 
a constant energy shift induced by this transformation. 

The dynamics of the system can be analyzed through 
the master equation, where the coherent electron trans- 
port through the molecules is associated with the off- 
diagonal operators s, s' . In the next section we explicitly 
give details of the master equation. 



III. MASTER EQUATION 

To obtain the master equation, we switch to the in- 
teraction picture. In this picture a^a^ and a R a R are 
unchanged since they commute with the non-interacting 
Hamiltonian: 



a\a L (t) = n L (t) = n L 
a R a R (t) = n R (t) = n R . 



(12) 
(13) 



For the term Ht, we define new operators that group the 
electron tunneling terms and vibrational mode operators: 



s (t) = pexp(iei)A t , 

fit (£) =pt exp (_ iei )X t t , 

where e = ei, — e R . Then 

H T (t) = T c [s (t) + St (t)] 



(14) 
(15) 

(16) 
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In the interaction picture, the equation of motion is 



= "ft (*)'/»(*)] 



Using these we write the following coupled integral equa- 
tions for the expectation values of Ul, n R , s, s\ 



[H T (t)+H v (t),p(t)}, 



(17) 



(n L ) t (h L ) = -j: J dt' [T c «S (t)) - <*t (*))) 



+r L (n L ) - T L f s } 



(24) 



and the density matrix can be obtained by iterations. By 
keeping terms to second order in the couplings to the the 
source and drain, and making the standard Born-Markov 
approximation, the density matrix p (t) can be approx- 
imated as p(t) ~ Pq° (t) <g> Pq os . Then, after tracing 
out the leads, the master equation can be written in the 
Lindblad form 



(n R ) t - (h R } 



dt' [T c ((S(t))-(5t (*)» 



% 

h jo 

-Tr (nil) - TrJd 



(25) 



dp DO (t) 
dt 



-^[H T (t),p 



DOl 



(s) t - (s) = jf df e «(*-*')T c 



-ffsV[a L ] 
^fnV[a R }, (18) 



{(n L X t xl) t -(n R xlx t )^ 

dt > e ie(t-t') /s( t >)xlx t \ 



(r L + iV '•' 



(26) 



where 



V 



ee a\a lP DO (t) 2a iP DO (t) a\ + p DO (t) a\a, _ (gt) Q = _ j f ^ e ~-(*-0 T( 



(19) 

27 [en] ee ^ojp 00 (t) - 2al^° (t) a, + p 150 (t) ^aj 

(20) 

and /s and /u are the Fermi distribution functions 
for the left and right reservoirs, respectively, fs = 



{(n R XtX v ) - (n L X t ,xf)) 



(r L + r R ) 



di J e - ie (t-t') /x t ,x}8* (*') 
(27) 



+ 1 



and f D = [e h ^-^) 



+ 1 



In all 



the later results we choose an infinite bias so that f s = 1, 
/d = 0. Since we are interested in the strong-coupling 
limit, one may also argue that the perturbative coupling 
to the leads should be derived after diagonalization of 
the combined electronic/vibrational system. This is an 
interesting avenue for future work. 



IV. EXPECTATION VALUES 

At this point the molecule-vibrational mode coupling 
is still described exactly. To achieve an analytical solu- 
tion it is convenient to rewrite the master equation as 
a closed set of equations for the expectation values of 
the various molecule operators, following the formulation 
used by Brandes^i to describe quantum dots coupled to 
a bath of oscillators. To do so we need the following 
commutators, 

[n L (t),H T (f)]=T c [s(t)-sHt]] (21) 

[s (t) , H T (*')] = T c e i£ (*-*') {n L xl,X t - n^X,} 

(22) 

[St (t) , H T (*')] = r c e ie (*-'') [n R X t ,x\ - n L X v x\} 

(23) 



This is still exact, but intractable. To make progress one 
assumes that the molecule operators and the oscillator 

ones are separable, so (n^XtXt) can be written as 



n L XtXl 



= (n L ) (1 - 
ee (n L ) F, 



e- nf}UB X t xl 



(28) 



with j3 — ksT. We have also assumed that the vibra- 
tional mode is in equilibrium due to contact with a ther- 
mal bath at temperature T. The other properties of the 
mode are ensconced in the two-time boson correlation 
functions T . In the appendix we give an analytical form 
for J- in the limit of an undamped single mode. Then, 
one can in principle solve the above equations of motion 
for arbitrary molecule-vibron coupling strength (which 
includes non-Mar kovian properties of this interaction). 

The polaron transformation method is powerful in that 
it allows us to gain simple analytical forms in some cases, 
and in general "plug-in" arbitrary bath correlation func- 
tions for J- although we need to be mindful of the regime 
where the method is valid. Some works specify that it 
gives qualitatively correct results for small T c , while, a 
comparison to the NIBA 58 suggests it is very accurate 
for e = or e > T Cl g. For the case of a single-mode, 
which we consider here, understanding when the results 
are accurate can be challenging, as we will discuss in the 
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next section. Generally, we do find them to be qualita- 
tively correct, though additional back-action effects arise 
in the exact numerical treatments. 

Finally, before showing explicit results, we make an 
additional assumption. The above equations allow for 
double electron occupation of the dimer molecule. Allow- 
ing such occupation is interesting, but prevents us from 
obtaining analytical results. Thus, as mentioned earlier, 
we impose another condition equivalent to the notion of 
Coulomb blockade in quantum dots. This is valid if the 
molecule is small and the bias not too large (i.e., this 
Coulomb repulsion is the largest energy scale). In this 
case only one electron can inhabit the dimer at a time 
and sequential tunneling occurs. 

Rather than explicitly include this effect with a non- 
linear repulsion against double-site occupation in the 
Hamiltonian, we can simply replace the fermionic opera- 
tors with projectors onto the single electron states 

n L = \L)(L\, n R = \R)(R\, (29) 

with 

n L +n R + \0)(0\ = l, (30) 

where |0)(0| is the no-electron state. Similarly, s = 
\L)(R\, = \R)(L\. The combinations of operators 
in the Lindblad terms describing the coupling to the 
leads and Equations [241127] are subtly altered by this 
change, to preserve the reduced Hilbert space. We do 
not give explicit details here as the derivation is equiv- 
alent to the previously studied model of a double quan- 
tum dot coupled to a phonon bat h 61 ' 71 , except where 
the single-mode-correlation function is replaced by the 
spectral function for a particular bath model. For com- 
pleteness we give the altered equations of motion in the 
appendix. It is these equations of motion, under this ad- 
ditional assumption, which give us the results in the next 
section. 



V. RESULTS 
A. Current 

The expectation value of the current from the source 
to the left molecule is II (t) = —eT^ (1 — fs) (nL), 
from the right molecule to the drain is I R (t) — 
— eT R (l — fo) (n R ), and from the left molecule to the 
right molecule is I L (t) = -eT c ((s (t)) - (st (t))). In the 
steady-state these three currents are the same, and can 
easily be evaluated from the expectation values of the 
molecule operators descried above. 

The most straightforward method to solve such equa- 
tions is to explicitly Laplace transform the equations of 
motion first and then find the coefficient of the Laplace 
parameter 1/z as z — > in the Laurent expansion of the 
expectation values^. We can also use the full z depen- 
dent forms to find the current-noise spectrum, as outlined 



in the next subsection. As an example, we explicitly solve 
here for the current using the J- function described in the 
appendix. As discussed earlier, this allows us to inves- 
tigate the unusual limit of strong coupling of the trans- 
port process to an undamped single mode. Note, in all 
of the following results, we set K = 1, e = 1 and typically 
Mm = 1, apart from in the noise power results. 

The top figure in Fig. [T] shows the current as a func- 
tion of the energy difference between the left and right 
molecules (scaled with the resonant frequency of the vi- 
brational mode) for various vibrational mode coupling 
strengths (z = 0.1 black- and z = 0.5 blue-dashed lines) 
and without the mode coupling (solid red line) with a 
thermal occupation of the vibrational mode of TV = 0.05. 
The bottom figure in Fig. Q] shows the current for dif- 
ferent temperature (N — 0.01 black- and N = 1 blue- 
dashed lines) while keeping the coupling strength con- 
stant at z = 0.3. 

When e is much larger than the tunnel coupling, the 
electrons cannot tunnel between left and right molecules 
and the current approaches zero. With the vibrational 
mode coupling "on" , the current peaks when the en- 
ergy bias e equals the vibrational mode resonant fre- 
quency and its multiples (as also observed in the case 
of quantum dots coupled to a single phonon modei 3 -). 
Positive e spikes correspond to phonon-emission-assisted 
transport and negative e spikes correspond to phonon- 
absorption. In figure[5]we see that increasing the coupling 
strength narrows the lower emission peaks and broadens 
the higher peaks. The bottom figure in Fig. [5] shows how 
raising the temperature broadens and raises the absorp- 
tion peaks. In addition, between peaks there is a strong 
suppression of the current, akin to the Franck-Condon 
blockade. If we directly observe, around the first reso- 
nance, the steady state occupations of the various states 
we find that this blockade effect strongly localizes the 
electrons in the left molecule. 

The amplitudes of the absorption peaks strongly de- 
pend on temperature. Expanding the expression for the 
current around the first absorption resonance we find 

-(V 1 i Tr i 1 + 2N V 1 (si) 

max y L 4T 2 f Tr ^ + ^ J ■ \ a > 

In contrast, the height of the emission peaks is only 
weakly dependent on the thermal occupation, and, ex- 
actly at the resonant point, the maximum is not depen- 
dant on the coupling strength or the order, e = nui, of the 
resonance. For low temperature the height of resonance 
is proportional to the zero-bias current 

/^L=^=o(l + / e =o), (32) 
4T 2 

/e =° = r R + 4T?(2r-i + r- L iy (33) 

However the width around the resonance is strongly de- 
pendent on both the coupling strength and the order 
nu of the resonance peak. In addition, the properties 
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FIG. 1: (Color online) (Top) The current In/e, as given by 
the polaron transformed results, through the molecule as a 
function of molecular energy bias e/tj m = {f-L — £R)/oJm- The 
solid line is the uncoupled case, while the dashed lines are 
for different magnitudes of the coupling to the vibrational 
mode: z = g/u) m , (z — 0.1 black- and z = 0.5 blue-dashed 
lines, with a thermal occupation of the vibrational mode of 
N — 0.05). Under this polaron transformation assumption, 
and for a completely undamped vibrational mode, we see ex- 
tremely strong vibrational emission resonances in the current. 
(Bottom) The current In/e for no coupling (solid line) and 
for coupling z — g/uj m = 0.3, and for several different choices 
of the thermal occupation of the vibrational mode N = 0.01 
black- and N = 1 blue-dashed lines. As we raise the thermal 
temperature of the mode, absorbtion resonances occur, and 
the linewidth of the emission resonances becomes broadened. 
In both figures we set Tl = I\r = 0.1, T c — 0.1, and ui m = 1. 



of each resonance peak are almost entirely defined by 
the appropriate term retained in the sum in Eq. (1501) . 
Thus, overall, we can say that in-between resonances 
there is a strong suppression of the current due to the 
strong coupling to the vibrational mode, and exactly on 
resonance (for emission) the transport channel is com- 
pletely transparent. This represents a unique kind of 
vibrationally-mediated transistor, where the transport 
can be exquisitely sensitive to small changes in the ener- 
gies of the electronic states in the dimer. 

The fact that the height of the peaks is independent of 



the coupling strength is surprising. Under exactly what 
regime are these results valid? Recall that we made two 
approximations in the derivation; we assumed that the 
electronic and bosonic operators were separable, and that 
the bosonic correlation functions were entirely described 
by their thermal equilibrium state, so that 



n L x t xy} = 



n L ) (1 - e 



-fat 



X t XJ 



(34) 

This implies a neglect of the back-action onto the single- 
mode dynamics. In addition, we assumed that the dy- 
namics of the bosonic system were not damped (though 
this can be included via numerical evaluations of the cor- 
relation function). These two approximations can be in- 
terpreted in the following way; that we have coupled the 
electron transport to a semi- classical oscillator that is 
capable of emitting and absorbing energy at very spe- 
cific frequencies. The fact that we neglect back-action 
does indicate there should be a quantitative difference to 
arise from the exact solution. In the final section we will 
compare these results to two exact numerical models. In 
addition, as mentioned earlier, the set of approximations 
made here is, at least in the many-mode case, formally 
equivalent to the non-interacting blip approximatio n 58 ' 59 . 



B. Noise spectrum 

The current-noise power spectral density is a useful 
and easy way, with this formalism, to gain some in- 
sight on the transient dynamics. The current noise 
power has three components: particle currents through 
left molecule Sl (w), the right molecule S r (uj) and the 
charge noise spectrum Scn 

S H = a 2 S L M + fi 2 S R (w) - a/3io 2 ScN (w) . (35) 

The coefficients a and (3, with a + /3 = 1, depend on the 
capacitances between junctions and molecule. Here we 
focus on the contribution from the right particle current, 
which can dominate if j3 3> a. This contribution can be 
obtained from the MacDonald formul a 72 ' 73 and 



S R (w) = 2oje 2 / dt sin {tot) — (n 2 (t)) - (t (I) 2 ) 
Jo at I V / . 

= 2eI{C(n t {t),h 3 {t + T)J 

+C(h i (t),n j (t + T))_. u }. (36) 

The steady state of Sr (oj) can be calculated as 

S R (w) = lei {1 + T R [h R (-tw) + h R (iw)]} (37) 

The frequency dependent n R (iu>) is again easily ex- 
tracted from the Laplace transform of Eqs. [M] and [25] , 
However, as noted in^i the validity of the noise spectrum 
formulation with the polaron transformation treatment is 
ambiguous can easily produce non-physical results. Fig- 
ure |31 show the spectral density with both no coupling 
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FIG. 2: (color online) (Top) The current I R /e through the 
molecule as a function of the molecular energy bias e/uj m 
and vibronic coupling strength z = g/ui m , with thermal oc- 
cupation N = 0.05. (Bottom) The current I R /e for coupling 
z = g/uim — 1 as a function of e and thermal occupation N. 
In both figures white is a large current and purple is zero cur- 
rent, and these results are from the polaron treatment, and 
again for 1^ = T R = 0.1, T c = 0.1, cu m = L 



(red line) and weak coupling for two different vibrational 
frequencies ui m (uj m = 1, black- and uj m = 2, blue- dashed 
lines) for T c = 0.2, r L = T R = 0.1, and N = 0.05. 
There is a clear large resonance from the bare Rabi os- 
cillations of the electron tunneling between molecules at 
uj = 2T C = 0.4, as well additional resonances at uj = u rn . 



VI. NUMERICAL COMPARISON: EXACT 
SOLUTION AND THE HIERARCHY EQUATION 

For small thermal occupation N, reasonable coupling 
strength z and by including damping of the vibrational 



1 05 



-£ 1.00 
if) 



| 0.95 

o 
a. 

a o 90 

v> 
o 
c 

■£ 0.85 

q> 
i 

O 0.80 



0.75 




05 



10 15 
Scaled Frequency, 



2.0 



2 5 



FIG. 3: (Color online) The current noise power Fano factor 
Fr(uj) = S R {uj)/2eI R , for F L = T R = 0.1, T c = 0.1, z = 
0, 0.3, and cj m = 1, 2, and iV = 0.05. There are two resonances 



w = 2T C and w 



as well a an additional higher resonance 



at 2w m . The latter two resonances are indicated by the two 
blue circles. 



mode (or an effective multi-mode environment) we can 
solve the molecule- vibration coupling exactly with a va- 
riety of different numerical approaches. Here we use the 
Hierarchy equation of motion for a Lorentz bath spectral 
density recently introduced by Ma et al§^. Essentially we 
take our original equation of motion, but write the inter- 
action with the vibrational mode as a bath of oscillators. 
We then assume that this bath has a Lorentz spectrum, 



1 



J(w) = - 



7T (OJ - UJ m ) 2 + 7 2 



(38) 



where A is equal to the square of the molecule- vibron cou- 
pling strength A = g 2 , ui m is the fundamental frequency 
of the vibronic excitation, and 7 is the broadening or 
dampin g of the vibronic mode. We then follow the steps 
in Ref. [6J] to reach the following hierarchical equation 
of motion (again setting h — 1), 

d 



k=l 

-4£ n * [V x + (-l) k V°] pn-sM (39) 
fe=i 

The superoperator notation introduced in this equation 
is V x . = [V, .] and V°. — { V, .}, where V is the electronic 
operator which describes the coupling to the Lorentz 
bath, which in this case is V = — tir. Here, ft is a two- 
dimensional index with positive integer elements identi- 
fying both the true system density matrix /9(o,o) and the 
auxiliary operators which encode the effect of the bath 
P(>o,>o)- The other vectors are defined as e% = (1,0), 
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e*2 = (0, 1) and fx = (7 — iuj ,-f + iuj ). £ is the super- 
operator describing the coherent and incoherent dynam- 
ics of the molecule transport. Note that again here we 
assume a large double-occupation repulsion (Coulomb 
blockade), so that the electronic and transport proper- 
ties are simply described by the single electron occupa- 
tion projectors, «l = \L){L\, ur = \R)(R\, s = \L){R\ 
and, 



0.03 



c\.\ = 



—I 



^(n L -n R ) + T c {s + «t) j 



2 D[S L] - r fD 



SI 



(40) 



where Sl = |0)(I/|, Sr — \0)(R\, and V is as defined in 

Eqs. dinHsnD- 

The tier of the hierarchy in this method must be trun- 
cated at some reasonable value, n = (N C ,N C ). Typically 
this is done by increasing N c until convergence is found. 
If A is increased or 7 is decreased the N c necessary for 
convergence rises (representing the increased occupation 
of the bath, and the increase in associated degrees of free- 
dom). The hierarchical equations of motion are solved by 
setting P(o,o) to the desired molecule initial state, setting 
all auxiliary density matrices to zero P(>o,>o) = 0, and 
employing direct numerical integration. Here we are pri- 
marily interested in steady-state results so we integrate 
until the system probabilities become stationary. The 
steady state can also be found by finding the null vector 
of the matrix defining these equations of motion. 

To check the validity of this method we also directly 
diagonalizc the equation of motion (without the polaron 
transform, and with an additional Lindblad damping of 
the vibrational mode with rate 7^ = 7/2) in the Fock ba- 
sis, and integrate the resulting coupled equations of mo- 
tion. We find that, given convergence, the results for the 
two numerical methods in the steady-state match exactly. 
Comparing to the polaron results we see some similarities 
and some differences. In the top figure of@]we see that 
the numerics, like the polaron treatment, predict equal 
height emission peaks for the different resonances as 7^ 
becomes small. However as 7^ is decreased further there 
is also a reduction in the overall current magnitude. This 
is an additional localization of the electron due to the 
strong coupling to the vibration not captured by the po- 
laron treatment. This phenomena is more clearly shown 
in the bottom figure of 2] where we show the current ver- 
sus the vibrational mode damping at the first resonance 
peak, e = uj m . We see that for all coupling strengths 
there is a common maxima following by a suppression 
of the current as 7;, is reduced. This indicates that, for 
the purpose of current enhancement, there is an optimal 
finite vibrational mode damping. 



VII. CONCLUSION 

In conclusion, we have analyzed the properties of co- 
herent electron transport through a dimer that is strongly 
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FIG. 4: (Top) The current Jjj/e as given by the exact numeri- 
cal treatment for z = 0.3. The red solid curve is for 75 = 0.005 
and the dashed blue curve is for 7^ = 0.0025. We find that as 
we reduce the vibrational mode damping current qualitatively 
behaves like the polaron treatment results. However this is 
an additional overall reduction in magnitude not predicted by 
the polaron results. (Bottom) The current Ia/e as a function 
of damping for e = ui m , for three different values of the cou- 
pling strength. We see that the three curves have a similar 
maxima, but are suppressed when 7;, is reduced (the z = 0.1 
curve is also suppressed for 7;, < 0.001). 



coupled to a vibrational mode. When the energy dif- 
ference between the electronic states in the two sites of 
the dimer is larger than the tunneling coupling, the cur- 
rent is suppressed. On the other hand, a semi-analytical 
polaron treatment suggests, when this energy difference 
is near the vibrational resonant frequency, an enhance- 
ment of the electron transport appears as extremely 
sharp and equal-height current resonances. This strong 
coupling also induces Frank-Condon blockade, suppress- 
ing the current between the resonances. The vibra- 
tional mode temperature determines the nature of the 
phonon-assisted transport. At low temperature, phonon- 
emission-assisted electron transport is the dominant fea- 



ture. As the temperature is increased, the phonon- 
absorption-assisted electron transport also starts to play 
a role. 

We also examined two exact numerical treatments of 
the same system. These two treatments were found to 
be identical in the steady-state, and qualitatively agreed 
with the polaron treatment. In addition, these results 
suggest that an additional localization or current sup- 
pression mechanism appears for very low vibrational 
mode damping, and thus there exists an optimal finite 
damping. 

Our investigation is primarily motivated by studies of 
electron transport through single molecules. However, 
our analysis is also applicable to artificial systems, in par- 
ticular quantum dots in carbon nanotubes, where strong 
coupling can occur. The connection to large systems like 
doubly-clamped semiconductor beams and cantilevers is 
more tenuous, as in those cases we expect the coupling 
to be weak and thermal phonon effects to dominate the 
transport properties. 
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X (t) and (t) operators become 



X (t) = cxp 
= exp 

x exp 



-z (tie 

(-z&V"*) exp (zbe- iut \ 
l r 

2 



= cxp —z\re %uit exp zbe 



tfe- tut ,zbe- %ult 

iujt _ — ■ 



(t) = e~^ 2 exp (zSV"*) exp (-zbe 



then, XtXj, becomes 



XtX], = e~ z e 



Putting the operators in normal order by using 



e iH T t/hQ e -iH T t/h _ q 
1 

+ 2! 



-tH T ,0 



:tH T , 



-tH T ,0 



+ 



(43) 



(44) 



(45) 



(46) 



and after noticing e~ M f (i) e M = e~ zb t' e zit e zb t' , with 



= e zbt , A = 



e- zb l> 



{be~ lut ^j e zi l' = e~ iut (b + ze lu)tl ) 



r 2 -w{t-t') 



(47) 



IX. APPENDICES 



Boson correlation functions 



so that 



e zht e zb l> = exp 



and we obtain 



r 2 -iw(t-f) 



z6+e ,0Jt ' ze ,0Jt b 

e e 



(48) 



Here we show how to derive the analytical form for 
the correlation function of a single undamped oscillator 
mode. This is equivalent to a recent derivation given in 
[o0l |. We start with a similarity transformation 



~ M f (i) 



e M = e- M 



£a„i) e M (41) 

n=0 J 

a n A = f (A 



E-M a M 
a n e Ae 



n=0 



n=0 



which suggests that one needs to consider only the trans- 
formation on each operator separately. Utilizing the in- 
teraction picture of the boson operators and their com- 
mutation 



-z&e iwt ,zbe~ 



(42) 



X t x\, = exp 



x exp 



-z 2 (l - e-M*-*') 
-zb (e~ iut ' -e~ iut 



exp 



(49) 



The correlation J- can be arranged in more a convenient 
form 



T= (l-e-^Y, 

n=0 



E 



\ > \n 

n=0 



(50) 



Let u — z (e~ iut ' - e"^'), u* = z (e 1 ^' - e iut \ and 
expanding the exponents in a power series and use the 
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property of destruction operators will condense the form the exponential part of Eq. (|52|) in terms of exp ±iwr 
to the Laguerre polynomial or order n 



{-I) 1 n 



L 



('"' 2 ) 

n(\u\ 2 ) (51) 



Using one of its generating functions 
E„=o L n (M 2 ) r = e-^H7(l - gives 



■F(r) = e- 



-z 2 (l+2iV) ^ ,2^ «+P ^ C 1 + N T c i(s-p) 

s,p=0 S 'P' 



(53) 

Physically, the s terms account for absorption and p for 
emission of phononsSs with frequency ui . Writing as 



J-(r) =exp{-z 2 [(l + iV)(l-e-^ r ) + N (l - e iuJT )]} 

(52) 

where r = t — t' . This form is still not easy to evaluate 

in the integral form of Eq. (|27[) . Therefore, we expand we can rewrite 



n=p - S = riemit - «ab, 
p = S + ™, 



(54) 
(55) 



J 7 (r) = exp —inujT exp — z 2 (1 + 27V) 

n= — oo 



1 + N\ 2 



N 



I n (2z 2 VWTiV)) 



(56) 



where I n is the modified Bcsscl function of the first kind. 
This form can be numerically evaluated. Similarly, the 
other correlation functions are 



(57) 
(58) 



B. Equations of motion for the Coulomb blockade 

limit 



Xl,X t ) = (x v Xt) = F (r), 
X\X t >)=T{r). 



(n R ) t - (n«) = - jf dt' [T c ((S(t)> - (it (t))) 

-I* <n R )] (60) 

(g) t -<6f) = -ijf t dfe"(*-*')r c 

dt'e 1 ^-*') (s(t')xl r X t \ (61) 



2 ./,, 



To obtain the analytical and numerical results we show 
in the second half of this work we imposed an additional 
condition on our model, that of strong Coulomb repul- 
sion against double occupation of the dimer. In this case 
we replace the operators with projectors onto a single 
electron basis, and obtain the following slightly altered 
equations of motion. 



(n L ) t - (n L ) Q = J o dt' [T c «S(t)) - (gt (*)» 
+r L (1 - (n L ) - (n R ))} 



(59) 



h Jo 



dt > e -Ht-f) T 



[(n R xlx v )-(n L X t ,x})} 



dt 



(62) 
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